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Characterised by long term cooling and abrupt ice sheet expansion on Antarctica —14 Ma ago, the mid 
Miocene marked the beginning of the modern ice-house world, yet there is still little consensus on its causes, 
in part because carbon cycle dynamics are not well constrained. In particular, changes in carbonate ion 
concentration ( [COa^"] ) in the ocean, the largest carbon reservoir of the ocean-land-atmosphere system, are 
poorly resolved. We use benthic foraminiferal B/Ca ratios to reconstruct relative changes in [COs^"] from 
the South Atlantic, East Pacific, and Southern Oceans. Our results suggest an increase of perhaps —40 fimol/ 
kg may have occurred between — 1 5 and 14 Ma in intermediate to deep waters in each basin. This long-term 
increase suggests elevated alkalinity input, perhaps from the Himalaya, rather than other shorter-term 
mechanisms such as ocean circulation or ecological changes, and may account for some of the proposed 
atmospheric CO2 decline before —14 Ma. 

The mid Miocene, including the mid Miocene climate transition (MMCT) at —13.8 Ma', represents one of 
the three most significant cooling episodes of the Cenozoic'^''', along with the Eocene- Oligocene transition'' 
(—34 Ma) and onset of Northern Hemisphere Glaciation (—2.7 Ma). Yet in contrast to these events, 
comparatively little is known about the causes and feedbacks of this mid Miocene cooling. Set within a long 
term cooling trend from 15 to 10 Ma, the rapid expansion of the East Antarctic ice sheet at 13.8 Ma''"^ is associated 
with a sea level fall in the order of —60 m*"'', and a coeval drop in deep ocean temperature of — 2°C^''*. This glacial 
expansion was accompanied by the extinction of Antarctic tundra', a — 6-7°C fall in sea surface temperature in 
the high latitude southwest Pacific'", the development of perennial sea ice in the Arctic" '^, the development of 
significant temperate terrestrial biotic provinces including the early expansion of C4 grasses and significant 
tectonic events that would have impacted on the global carbon cycle (see Supplementary Information). 
Hypotheses for the causes of this transition include CO2 drawdown, perhaps from a cessation of mid Miocene 
volcanism" or enhanced terrestrial weathering and erosion'^, priming the climate system for ice sheet growth 
during low amplitude orbital eccentricity at 13.8 Ma'. Alternatively the onset of a deep Antarctic Circumpolar 
Current may have thermally isolated Antarctica due to tectonic subsidence in the Scotia Sea"". 

There have been numerous efforts to reconstruct atmospheric CO2 through the Miocene''*"", and in contrast 
to previous suggestions that temperature may have become decoupled from CO2 during the mid Miocene"''*, a 
recent study suggests a decline of — 100 ppmmay have occurred between —15 and 14 Ma before the MMCT'"*. As 
surface ocean CO2 is in approximate steady state with the atmosphere, and the relatively small surface ocean 
reservoir is continuously replenished by the much larger deep ocean reservoir (>1 km), the speciation of 
dissolved inorganic carbon (DIC) of the deep ocean is linked to atmospheric CO2 (refs 20,21) and therefore a 
— 100 ppm decrease in CO2 should be documented within deep sea sediments. In the ocean, DIC occurs as several 
species that maintain acid/base equilibrium: 

DIC = [CO2 (aq)]^[H2C03]^ [HCO3" 1 ^ [CO3" 1 



As the ratio of ocean water alkalinity (ALK; the sum of bases in solution, in sea water predominantly: 
[HCO3"] + 2[C03^"] + [B(OH)4"]) to DIC increases (for example from increased continental weathering 
or increased biological productivity), the proportion of the DIC pool present as sea water [COs^"] increases, 
thus shifting the balance away from [C02(aq)] and potentially drawing down atmospheric CO2 to recover 
equilibrium^"'^'. Previous attempts to constrain past ocean carbonate chemistry have focused on reconstructions 
of the calcite compensation depth (CCD) linked to [C03^"], which are of low resolution and subject to 
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Figure 1 | Global palaeogeographic map showing the continental configuration at 14 Ma. The positions of ODP Sites discussed in this study are 
indicated as filled circles. Possible major sources of deep Northern Component Water (NCW), Southern Component Water (SCW), and Tethys-Indian 
Saline Water (TISW), are indicated by schematic arrows. MoUweid projection of modern continents (red) is shown on a palaeogeographic 
reconstruction, generated from^^, of continental plates (grey) centered at 14 Ma. 



significant complications"'^, or reconstructions of [COs^"] directly 
with benthic foraminiferal shell Li/Ca and Mg/Ca', also subject 
to complications such as a temperature effect' (Supplementary 
Information). 

Recent studies^''^"* on core top samples from several ocean basins 
have shown a strong quantitative relationship between benthic for- 
aminiferal B/Ca ratios and deep water [C03^"], which has been 
further supported by subsequent down-core studies^^ Here, we 
reconstruct deep and intermediate water [C03^"] from six Ocean 
Drilling Program (ODP) core sites in three ocean basins (Fig. 1) 
using B/Ca of the epifaunal benthic foraminifera Cibicidoides mun- 
dulus (Fig. 2, see Methods). The species C. mundulus is unique in 
being long-ranging, allowing the use of this proxy as far back as the 
mid Miocene. 

Results 

Cibicidoides mundulus B/Ca. In the ~1 Ma lead up to the MMCT 
(~15.5-14 Ma), benthic foraminiferal B/Ca ratios from all six sites 
increase (Fig. 2, Supplementary Table SI), indicating a probable 
increase in deep water [COs^"]. All samples have low Al/Ca values 
well below 100 |imol/mol (Supplementary Table SI), indicating little 
potential clay contamination''^''. In the South Atlantic (Fig. 2b) the 
deepest Site 1226 (—3.6 km palaeo-water depth) shows a larger 
increase in B/Ca over this interval than shallower Site 1264 
(2.3 km), indicating a possible greater increase in [COa^^] in deep 
waters relative to intermediate waters. In the East Pacific (Fig. 2d) 
and Southern Ocean (Fig. 2e), the magnitude of the B/Ca increase at 
both the shallow and deep sites appears to be comparable, with a 
slightly larger overall increase in the East Pacific perhaps due to a 
modified ocean circulation pattern. In the ~ 1 Ma after the MMCT 
(-14-13 Ma) deeper sites in the East Pacific (Site 1237, 2.8 km) and 
Southern Ocean (Site 1 168, 2.3 km) exhibit further increased benthic 
B/Ca, whilst at nearby shallower Sites 1236 (0.9 km) and 1171 
(2 km) benthic B/Ca is seen to drop slightly. This B/Ca decrease at 
the shallower sites may be explained by either a decrease in global 
intermediate water [C03^"] perhaps due to changes in ocean circula- 
tion, or more likely, localised depletions in [COa^^] due to expansion 
of the oxygen minimum zones associated with enhanced primary 
productivity (both sites are situated in proximity to the modern 
depth of [C03^"] and oxygen minima; Supplementary Fig. SI). 

Miocene oceanic [COs^"]. B/Ca ratios from all sites were converted 
to [C03^"] using the 0.69 |imol/mol per |j,mol/kg sensitivity 
obtained from a global core-top calibration^'', past water depth 
estimates, and mid Miocene estimates of seawater boron and 



calcium ratios (B/Ca^w) and concentrations (Supplementary 
Information). There is yet to be a consensus on mid Miocene 
boron and calcium concentrations""*'^' and future estimates may 
therefore change, although the relative changes in [COa^^] inferred 
from our records (Fig. 3) may not be greatly affected. It is likely that 
the B/Cas„ ratio changed over the interval 15.5-13 Ma, even 
considering the long residence times of boron (20 Ma)™ and 
calcium (1.1 Ma)" in the oceans. It is difficult to estimate these 
potential changes, but the dominant source of oceanic boron is 
considered to be continental discharge from weathering of silicates, 
carbonates and evaporites'^, and changing Himalayan erosion has 
been identified as an important variable on Cenozoic boron flux^". To 
test the sensitivity of our [COs"""] estimates to a changing flux of 
boron to the oceans, we consider the effect of a 10% increase in 
B/Cas„ given that sedimentation rates from Asian basins are 
estimated to have increased by —10% due to Himalayan erosion 
between 15.5 and 13 Ma'^' (Supplementary Fig. S2). A 10% 
increase in oceanic boron whilst calcium remained stable 
(conservatively assuming the majority of weathering products were 
silicates"), would have increased B/Cas„ by —3 |imol/mol, 
impacting [COs^"] estimates by only —4 |imol/kg (included in the 
error on Fig. 3). Calcium concentrations of the oceans between 15.5- 
13 Ma are actually estimated to have increased in the order of 50%^', 
largely from exposure of shelf carbonates as sea level fell particularly 
after —13.8 Ma^'*. Since the B/Ca of the exposed carbonates is 
unknown, the impact on B/Ca^w is unknown. However if calcium 
concentrations drove B/Ca changes, then we would expect to see 
similar changes in Mg/Ca and Li/Ca (Supplementary Table SI), 
which we do not. With the exception of the two data points from 
the shallowest sites (Sites 1236 and 1171) showing a decrease at 
— 13.5 Ma, we estimate deep water [COa^"] to have increased by 
an average of perhaps —40 |imol/kg between —15.5 and 13 Ma, 
with the largest increase occurring between —15 and 14 Ma (Fig. 3d). 

Over the interval —15.5-13 Ma, deep water CaCOs preservation 
may have increased in the South Atlantic, East Pacific and Southern 
Ocean (Fig. 2b, d and e), possibly reflecting the increasing [C03^"] 
(from B/Ca) and inferred deepening of the CCD. Carbonate preser- 
vation proxies in the North Atlantic (Fig. 2c, wt% CaC03), Indian 
Ocean (Fig. 2f, % CaCOs coarse fraction). West Pacific (Supple- 
mentary Fig. S3, CaCOs fragmentation) and Southwest Pacific 
(Supplementary Fig. S3, wt% CaCOs) show a long term increase, 
suggesting possible increases in [COs^"] in these locations also. 
There are significant uncertainties with using these records to infer 
preservation, as we do not yet have detailed sedimentation rate or 
CaCOs rain rate data^^, but taken together the CaCOa trends are 
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Figure 2 | Summary of ocean carbonate proxies. Deep-ocean B/Ca records 
(coloured symbols) compared with deep-ocean sedimentary CaCOs data 
(grey symbols) and deep-ocean oxygen isotopes (blue line) over the interval 
— 15 to 13 Ma. Open symbols indicate probable interglacial samples. 
Increasing values of B/Ca and CaC03 data indicates elevated bottom water 
[COs^^]. (a) Deep ocean S'^O from Site 1237\ indicating long term cooling 
with the sharpest drop at —13.8 Ma during the mid Miocene climate 
transition, (b) B/Ca with CaCOs wt% data from the Walvis Ridge^'. 
(c) CaCOa wt% data from the northern Mid-Atlantic Ridge", (d) B/Ca 
with CaCOs accumulation rate data from the Nazca Ridge^'"^^. (e) B/Ca with 
CaCOs wt% data from the Tasmanian Margin and South Tasman Rise*, 
(f) CaCO, coarse fraction data from the Wombat Plateau^. 



fairly consistent. The consistency of our [C03^^] reconstruction and 
CaC03 preservation records (Figs 2,3d) appears to point to a global 
increase in [CO3''"] of deep (>1 km) waters particularly from —15- 
14 Ma, as our records cover a large depth range (—0.9-3.8 km) and 

Age (Ma) 
14 14.5 




14 14.5 
Age (Ma) 

Figure 3 | Composite seawater [COs^"] estimates against other mid 
Miocene climate records, (a) Atmospheric CO2 reconstructed from 
planktonic foraminiferal 8"B (fiUed circles)'", and fossil leave stomatal 
frequency (open circles) (b) Deep ocean benthic foraminiferal 8'*0 from 
ODP Site 1237', indicating long term global cooling and ice-sheet 
expansion during the mid Miocene, and rapid cooling and ice sheet 
expansion at —13.8 Ma. (c) Deep ocean 5''C from Site 1237'. (d) Deep- 
ocean [COj^^] estimated from foraminiferal B/Ca data of several sites 
(Fig. 2), using estimates for mid Miocene B/Cas„ values and calculated 
palaeo-water depths (see Methods). Error largely associated with 
calibration uncertainty^' and possible changes to B/Ca^w due to continental 
weathering (see Methods). Blue symbols represent glacials, red 
interglacials, and green intermediate. Blue line is a best fit 5-point 
smoothing spline, with ± Is.d. of dataset (light blue). Circle indicates two 
samples that may have been affected by local oxygen minimum zone 
(OMZ) expansion, (e) Chemical weathering index from ODP Site 1148, 
South China Sea, as the ratio of chlorite/(chlorite + haematite -t- goethite) 
(Crat)'''- Lower values may represent increasing monsoon intensity over 
Southern China and associated intense weathering of the Himalaya". 
Other proxies (see Supplementary Information) indicate increased 
Himalayan weather and erosion during the mid Miocene. 
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geographic distribution of the major ocean basins (Fig. 1), and 
include both glacial and interglacial measurements (Fig. 3d). A study 
from ODP Site 761 at 2.2 km water depth in the Indian Ocean also 
concluded that there was an increase in local [C03^"] between —15 
and 14 Ma'. 

We recognise that there is substantial variability within our data- 
set, partly because the measurements come from various sites and 
from various times corresponding to both glacial and interglacial 
intervals and oceanic [CO3"'"] may be expected to have decreased 
during interglacials^'. However the long term increase in [COs^"] 
from —15 to 14 Ma presented here is superimposed onto these 
shorter time-scale variations. Significant basin partitioning of 
CaCOs may not be expected since similar values of benthic §"C in 
both the Pacific and Atlantic Oceans between 15 and 12 Ma 
(Supplementary Fig. S4) suggests a relatively homogenous deep 
water mass may have been present in the Southern Ocean"* and 
possibly the South Atlantic. Modelling experiments suggest that 
ocean circulation was dominated by Southern Ocean deep water 
mass formation^^, and sedimentological proxies document increas- 
ing influence of a deep Pacific southern water mass after —15 Ma"". 
Miocene Tethys Indian Saline Water, a significant component of 
deep and intermediate water bathing the Indian and Pacific 
Oceans (Fig. 1) and possibly warming Antarctica, appears to have 
diminished by —15 Ma"*' likely due to the shoaling of the Tethys 
Ocean basin'''. Northern Component Water is not thought to have 
played a prominent role in ocean circulation until after — 12 Ma''^'^", 
possibly due to the shoaling of the Greenland-Scotland Ridge. 
Neodymium isotopes from Walvis Ridge (used as an ocean circula- 
tion tracer) show no increased stratification and therefore North 
Atlantic deep water formation before —10.7 Ma". 

Discussion 

A long-term, whole ocean increase in [CO,^"] at > 1 km water depth 
and apparently enhanced CaC03 preservation (Fig. 2), particularly 
apparent between —15 and 14 Ma (Fig. 3d), implies an increase in 
the whole ocean ratio of alkalinity to dissolved inorganic carbon 
(ALK : DIG), as a charge imbalance driven by ALK is balanced by 
the speciation of DIG away from [GO2] towards [GOs^"] (ref 21). 
ALK is largely supplied to the ocean by continental weathering of 
emergent carbonates and silicates on land, and removed from the 
ocean by burial of GaGOs on the sea floor^". DIG is added and 
removed by these processes, but also removed by burial of organic 
carbon (Goig) in ocean sediments, and by air-sea exchange. There are 
several ways of causing increases in deep ocean [GOs^"], including a 
lowering of the GaGO^/Gorg rain ratio of primary export produc- 
tion^", increased shelf to basin fractionation from sea level fall and 
exposure and erosion of carbonate shelves*-''"'^'''^''''", changes in ocean 
circulation^"'^'', increased high latitude nutrient utilization increasing 
the [G03^"] of preformed surface water^"'", and increased contin- 
ental weathering'^"'"'"'*^. Processes such as ocean circulation and 
nutrient utilization can produce large short-term excursions in sea- 
water [G03^"], but are less likely to account for the long-term change 
observedfrom — 15.5 to 13 Ma. This is because carbonate compensa- 
tion, the process by which an elevation in deep water [GOs^"] 
increases the area of GaG03 burial on the sea floor and the sub- 
sequent removal of ocean water [G03^"], would bring the ocean to 
the initial steady-state on a timescale of several thousand years*^. 

Increased shelf to basin fractionation was proposed as the major 
cause of the widespread GGD fall during ice sheet expansion at the 
Eocene-Oligocene transition, as sea level fall drove carbonate depos- 
ition deeper into basins and raised oceanic §"0*'''°. However there is 
no evidence for a long-term sea level falP or long term 8"C increase 
(Fig. 3c) prior to 14 Ma, when we infer a notable [GOs^^] increase 
(Fig. 3d). If a large, long-term decrease in the GaG03/Gorg ratio were 
to account for some of the [GOs^"] variability from —15.5-13 Ma, as 
has been suggested to account for the short-term 8"C maxima at 



— 13.7 Ma**, we would expect to see increasingly organic-rich sedi- 
ments in marine sediment cores. In fact, the interval —16-14 Ma is 
associated with episodic increases in organic rich sediments around 
the Pacific margin*^ that are associated with global 5"G and S'^O 
maxima (Fig. 3), and higher productivity along the west coast of 
Africa""" before —14 Ma, indicating that changes in the GaGOs/ 
Gorg ratio may account for part the observed changes in [GOs^"] 
particularly before 14 Ma. However, no long-term secular increase 
in 5"G, consistent with a sizeable increase in productivity, is appar- 
ent in any records (Fig. 3c, Supplementary Fig. S4). 

We believe that the most plausible explanation for the long-term 
increase in [GOs^^] suggested by our records, particularly before 
14 Ma, was increased continental weathering and erosion'^. Uplift 
and erosion of the Himalaya, which has been linked to Genozoic 
climate change via the erosion and weathering of silicates which 
consumes GO2 (refs 33,34), would also have added ALK to the ocean. 
The Greater Himalaya experienced enhanced exhumation rates dur- 
ing the mid Miocene"" (Supplementary Fig. S2c), and plate config- 
uration reconstructions indicate the Himalaya met with the 
Intertropical Convergence Zone of high precipitation in the mid 
Miocene*"'. Several mineralogical proxies associated with monsoonal 
intensity also point to an increase in Himalayan regional chemical 
weathering in the mid Miocene'^ (Fig. 3e). Mounting evidence in 
support for enhanced weathering rates, and ultimately delivery of 
ALK to the oceans, is consistent with a peak in bulk sediment accu- 
mulation rates in the major clastic basins around Southeast Asia'^ 
(supplied by the Himalaya), and high fluxes of Sr from the Ganges- 
Brahmaputra Rivers*' during the mid Miocene (Supplementary Fig. 
S2e). Although dating is not yet able to precisely constrain the tim- 
ings of all these events, if enhanced silicate weathering and erosional 
deposition from the Himalaya was responsible for global cooling 
during the mid Miocene, as has been previously suggested'^'''''*^, then 
an increase in ocean ALK and hence [GOs^"] would be expected to be 
coeval. Other possible terrestrial sources of ALK were the Tethys 
Ocean, which underwent uplift and intermittent closure in the mid 
Miocene**, and the East Africa Plateau which also experienced peak 
uplift in the mid Miocene*' (Supplementary Fig. S2c), although the 
timing and magnitude of these events are even less well constrained. 
Whichever the source, enhanced continental weathering, consuming 
atmospheric GO2 and transporting ALK to the ocean, is consistent 
with our interpretation of elevated seawater [GOs^"] during the mid 
Miocene. 

In summary, our new reconstruction of deep ocean [GOs^"] 
(>1 km water depth) provides a tantalising insight into potential 
changes in the carbonate system during the mid Miocene. 
Comparison of these records with published 5"*0, 5"G, and CO2 
proxy records (Fig. 3, Supplementary Fig. S5) indicates enhanced 
input of ALK to the oceans as a possible contributor to mid 
Miocene climate change before the MMGT at 13.8 Ma. A likely 
source could have been the Himalaya which experienced increased 
uplift, weathering and erosion centred at —15-14 Ma (Supple- 
mentary Fig. S2). Whilst our records demonstrate that carbonate 
chemistry in the world's oceans changed during the mid Miocene, 
they also highlight the need to better resolve deep and intermediate 
water [COs^"] to fully understand the roles of weathering, ocean 
circulation and the carbon cycle during this complex episode of 
climatic change. However, taking our data at face value, we per- 
formed simple sensitivity calculations of the effect a hypothetical 
addition of ALK to achieve the estimated [G03^"] increase of 
40 |.imol/kg (Fig. 3d) may have had on ocean water CO2, using 
modelled estimates of initial Miocene ocean water ALK, DIG and 
temperature (Supplementary Information). Although necessarily 
subject to significant uncertainties such that absolute values must 
be considered speculative, if the deep ocean [COa^"] increase was 
mixed to surface water, it may have drawn down atmospheric GO2 in 
the region of — 16-25% (Supplementary Table S2). This is broadly in 
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the range of atmospheric CO2 reconstructions from fossil leaf sto- 
matal frequency^^ and planktonic foraminiferal 6^^B (Fig. 3a)^^. 
Further, in agreement with 5^^B reconstructions, the timing of our 
increase in deep ocean [C03^~] and putative drawdown in atmo- 
spheric CO2 is consistent with the largest drawdown in CO2 taking 
place between —15 and 14 Ma before the abrupt MMCT at 
13.8 Ma^^, and comparison of these two records does indicate a 
possible coupling between atmospheric CO2 and deep ocean 
[C03^~] (Supplementary Fig. S5). Atmospheric CO2 may have also 
been affected by other mechanisms such as volcanism^^'*^, and our 
reconstructions contain too much uncertainty to indicate how much 
CO2 was removed by the oceans. However, the general increase in 
deep ocean [C03^~] over the mid Miocene presented here does indi- 
cate a drawdown of atmospheric CO2 into the ocean may have 
occurred, and highlights the importance of constraining deep ocean 
changes when assessing mid Miocene carbon cycle dynamics. 

Methods 

Analytical procedures. Deep sea sediment samples of —40 cc volume were washed 
through a 63 |j.m sieve with deionised water, and oven dried at <30°C. A total of 10- 
25 specimens of the benthic foraminifera Cibicidoides mundulus were picked from the 
250-350 \xm fraction in each sample, and cleaned using the standard 'oxidative' 
treatment protocoP*''^^ Elemental/Ca ratios were analyzed by inductively coupled 
plasma mass spectrometer (ICP-MS) according to methods described elsewhere^^-^\ 
For the ICP-MS, the B blank is <2% of the consistency standard (B/Ca - 150 \imo\l 
mol). Recent detailed work on the foraminiferal genus Cibicidoides indicates that 
different morphospecies produce a certain degree of variation in B/Ca (ref. 24). To 
minimise possible error, we strictly selected only 'typical' C. mundulus using a 
conservative species concept (see Supplementary Information), from the 250- 
350 [im size fraction, and a relatively high number of specimens (>20 individuals 
where available). Replicate analysis carried out on two samples showed variability of 
<4 [j.mol/mol, in line with published replicate analyses of standards and samples 
giving precisions of <2% (RSD) for B/Ca^^. 

Miocene [COs^"] estimates. Values for A[C03^~] were calculated by modifying the 
equation: B/Ca = aAfCOj^'] + b (ref. 23), where a = 0.69 and /? - 119.1. To estimate 
mid Miocene B/Cas^ values (modern B — 415 |imol/kg; modern Ca — 0.01 mol/kg), 
we assumed B to be —1.07 times the modern value^^, and Ca to be —1.12 times the 
modern value^^. The resulting value for Miocene B/Ca^y^ is —96% of modern, which 
was used to modify the equation: B/Ca(Miocenc) = 0.96 * aAlCO^^'] + 0.96 * b. The 
calculated A[C03^"] was converted to [COj^"] by adding the [C03^"]sat of the 
pressure equivalent to estimated water depths at —14 Ma. See Supplementary 
Information for further details. 

Assessment of uncertainty. The error bars associated with B/Ca values (Fig. 2) 
include an analytical precision of ± 1.5 |j,mol/mol, and a replicate error of ±2 |j.mol/ 
mol. The error bars associated with [C03^~]sw (Fig- 3d) include the analytical 
precision of ±2.2 |j.mol/kg, the replicate error of ±2.9 |imol/kg, a calibration 
uncertainty of ±10 |imol/kg (2 a)^^, a sea level uncertainty of ±0.5 |j.mol/kg (from 
—50 m change), and a mid Miocene AB/Cag^^ uncertainty (from increased terrestrial 
erosion) of ±2 |j.mol/kg. 
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